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Abstract 

A pseudo Nambu-Goldstone boson as curvaton avoids the r^-problem of inflation 
wliicli plagues most curvaton candidates. We point out that a concrete realization of 
the curvaton mechanism with a pseudo Nambu-Goldstone boson can be found in the 
supersymmetric Peccei-Quinn mechanism resolving the strong CP problem. In the 
flaton models of Peccei-Quinn symmetry breaking, the angular degree of freedom 
associated with the QCD axion can naturally be a fiat direction during inflation and 
provides successful curvature perturbations. In this scheme, the preferred values of 
the axion scale and the Hubble parameter during inflation turn out to be about 10^'^ 
GeV and 10^^ GeV, respectively. Moreover, it is found that a significant isocurvature 
component, (anti) correlated to the overall curvature perturbation, can be generated, 
which is a smoking-gun for the curvaton scenario. Finally, non-Gaussianity in the 
perturbation spectrum at potentially observable level is also possible. 



1 Introduction 



The primordial curvature perturbation is caused presumably by some scalar field, which acquires 
its perturbation during inflation. For a long time it was generally agreed that the curvature- 
generating field would be the infiaton PJ. Then it was suggested instead that this field is 
some other 'curvaton' field [21 0] (see also which generates the curvature perturbation 
only when its density becomes a significant fraction of the total. The curvaton proposal has 
received enormous attention because it opens up new possibilities both for model-building and 
for observation.^ 

In all cases, the field responsible for the curvature perturbation must be light during infla- 
tion, in the sense that its effective mass rries is much less than the Hubble parameter H. (To 
be precise, we need during inflation m^^ < O.IH^, so that the spectral tilt 1 — n ~ ^mlg/H'^ 
satisfies the observational constraint.) If the responsible field is a curvaton, it should preferably 
also remain light after inflation, until H falls below the true mass Y]. These requirements are 
in mild conflict with the generic expectation from supergravity, that the mass-squared of each 
scalar field will be at least of order [H^ . This is the famous r^-problem of inflation. Ways have 
been proposed to keep the responsible field sufficiently light P^, the most straightforward of 
them being to make it a pseudo Nambu-Goldstone boson (PNGB) (for the curvaton see |21E])- 

For economy, and also to facilitate contact with observation, one would like a candidate 
for the responsible field to be one that is present in an already-existing model, designed for 
some purpose other than the generation of the curvature perturbation. Several such candidates 
have been proposed jHl CHI HH Il2l HH] , but in general they do not come with a satisfactory 
mechanism for keeping the curvaton sufficiently light. The purpose of this paper is to suggest 
a natural curvaton candidate, which is present in flaton models of Peccei-Quinn symmetry 
breaking, and which is a PNGB. 

The Peccei-Quinn (PQ) symmetry provides a nice solution to the strong CP problem It 
is an anomalous global symmetry, U{1)pq, spontaneously broken at an intermediate scale fpq 
with preferred value around 10^^ GeV providing enough dark matter (axion) of the universe. Its 
PNGB is the axion, which must be extremely light to satisfy the observational constraints J3] . 
In the context of supersymmetry (SUSY), two complex fields (or more) are needed to implement 
a global symmetry due to holomorphicity of the superpotential (unless the symmetry in question 
is R-symmetry). Hence, in SUSY, to spontaneously break the PQ symmetry one needs at least 
two complex fields, corresponding to two radial and two real angular degrees of freedom. One 

^More recently still it has been suggested that the curvaton acts by causing inhomogeneous reheating or 
through a preheating mechanism 
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of the latter is the axion, and the other is our curvaton candidate.^ We can define a symmetry 
(exphcitly broken of course) acting only on the combination of phases which corresponds to the 
curvaton, and then the curvaton is the PNGB of that symmetry. 

To understand how our model works, recall first that two fundamentally different paradigms 
exist for the PQ symmetry breaking fields. In one of them, the potentials involve only renormal- 
izable terms in the superpotential plus small soft SUSY breaking terms, leading to a normal 
Mexican-hat potential whose height and width have the same scale fpq. According to this 
paradigm, PQ symmetry breaking persists in the limit of unbroken SUSY so that the axion has 
a well-defined supersymmetric scalar partner called the saxion (as well as a fermionic partner 
called the axino). The saxion is the only light degree of freedom apart from the axion, its mass 
coming from SUSY breaking and being of order TeV for gravity-mediated SUSY breaking. The 
other two degrees of freedom are heavy with mass of order /pg. 

According to the other paradigm, the PQ symmetry breaking fields are instead fiaton fields 

[T6] . so-called because they are symmetry-breaking fields, which correspond to fiat directions of 

the potential (i.e., directions in which the quartic term is negligible). Their potential contains 

only soft SUSY breaking terms and non-renormalizable terms, which means that it is very fiat. 

In the limit of unbroken SUSY there would be no spontaneous breaking of the PQ symmetry. 

A nice feature of this paradigm is that the intermediate axion scale fpQ is not put by hand but 

is generated through a geometric mean of the SUSY breaking scale and the Planck scale. The 

other degrees of freedom accompanying the axion get their mass only from SUSY breaking, 

making them of order TeV for gravity-mediated SUSY breaking. In this paper, we point 

out that this kind of model contains a natural curvaton candidate as the angular degree of 

freedom. The radial degrees of freedom are less suitable for this purpose, because in the early 

Universe they presumably acquire the mass of order H which is generic |HI for scalar fields in a 

supergravity theory. In contrast, a mass of order H is unlikely to be generated for the angular 

degree of freedom, because it would correspond to a generalized A-term which is forbidden if 

the fields responsible for the energy density are charged under certain symmetries jHj. We will 

analyze how the model parameters like the axion scale, the Hubble parameter and the curvaton 

decay temperature are constrained to produce successful curvature perturbations. 

^By 'degree of freedom' we mean as usual a normal mode of the coupled oscillations of the fields, which after 
quantization corresponds to a particle species. 
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2 Supersymmetric realization of the PQ and curvaton 
mechanisms 



The minimal model of the type that we are considering is a variation of the non-renormalizable 
superpotential [El UHl HH] , which introduces two fields P and Q charged under PQ symmetry 
and one singlet S allowing the following superpotential 

W = h^-—HiH2 + \i^, + A2^ (1) 

where the U{1)pq and Z]^ charges are assigned as 

H2 P Q S 

U{1)pq: \{n + l) \{n + l) -1 +1 (2) 
: 1 1 \ a 

with a = e^'^*/^ and N = n + ?>. The field S is supposed to have a negative soft mass-squared 
of order the gravitino mass m^j2- This forces all the fields P, Q and S to develop the vacuum 
expectation values of the order 

/ \l/n+l , , 

~ (m3/2M^) . (3) 

This sets the intermediate PQ scale 

/pQ ~ 3 X 10^°, 1 X 10l^ 3 X 10^^ GeV (4) 

for n = 1,2,3, respectively, and ^3/2 = 300 GeV. The negative mass-squared of the S field 
can come from the renormalization effect [T^ or from the initial condition of soft masses. We 
do not assume the negative mass-squared for the P and Q fields, which would drive the scalar 
potential unbounded from below when some of the fields P, Q or S are set to zero. 

Note that the superpotential of Eq. (pj) provides a natural solution to the fi problem jSO] 
as the first non-renormalzable term leads to the right order of magnitude for the Higgs mass 
parameter; 

The above interaction will be the main source of the decays of the particles corresponding to 
P, Q or S into the standard (s)particles. 

The model under consideration contains six gauge singlet scalar particles. Three radial fields 
get masses of order m3/2 from soft supersymmetry breaking. Among three phase fields, ifp, 
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ifQ and ^95, one combination, yjp + yjg — 2(^s^ gets the mass of order m3/2 as the radial fields. 
The combination, (fp — v^q, is the axion. The other combination is our curvaton candidate, 
ifp + ifQ + ips, which becomes massive due to the soft supersymmetry breaking A term; 

Vsoft = AXi—^ h AX2--r- + h.c. , 

inducing also in the vacuum a curvaton mass of order 1713/2- (For a detailed calculation in 
a similar model of the mass spectrum, see [21J.) Note that, in our scheme, the curvaton 
parameters like its amplitude and decay temperature are fixed by the PQ scale and the /i term, 
as will be shown explicitly. 



3 The potential in the early Universe 

In order to consider the basic features of the curvaton field a associated with the QCD axion, 
let us simplify the original superpotential term in Eq. as follows; 

1^(0) = —^^, (5) 

where the complex field may be thought of as containing the phase field a and one radial 
component. This would be literally correct if the vacuum expectation values (VEVs) of P, Q 
and S happened to be equal, allowing the parameterization, 

= |0|e*^ = 101 exp{icT/V2v) (6) 

with V = (0) and a = \plvQ. (As noted in the Conclusion, our curvaton model can also be 
implemented using only 0, though one then loses the connection with the QCD axion.) 

With the above superpotential the general scalar potential in the early universe can be 
written as 

j0|2"+4 



F=(3C^i/^-m^)|0|^ + 



{CaH^A) +h.c. 



+ ^'^^ ' (7) 



where and A are soft supersymmetry breaking masses at zero temperature. Note that we 
have put negative mass-squared for the field at zero temperature to generate the PQ scale as 
discussed previously. 

The terms with H are expected to arise from the supersymmetry breaking effect during and 
after infiation jH], except during radiation domination [22] ■ The generic expectation for C,^ is 
\C^\ ~ 1, and we are going to assume C,/, ~ — 1 so that there is PQ symmetry breaking. As 
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is well known, a value \Ca\ ^ 1 can be achieved provided that relevant field values are small 
on the Planck scale, if a symmetry forbids low-order terms in the Kahler potential which are 
linear in the field responsible for the energy density [H]. To be safe though, we need an estimate 
of the actual value of Ca, considering for completeness also C^. 

During infiation, the supergravity potential is taken to be generated by a single F-term F/; 



V 



\Fi\ 



Ml 



(8) 



where 

= (^diW + WdiK/Ml) K^^* (diW + WdiK/Ml) * (9) 

and K^^* is an element of the inverse of the Kahler metric Kjj^ = didj^^K with J running over 
all of the scalar fields. Depending on the model, the field I might be the infiaton, the waterfall 
field or some other field [T]. 

We make the following assumptions; (i) in the Kahler potential, terms linear in I are negli- 
gible (this can be a consequence of a global or gauge symmetry), (ii) the field / is much smaller 
than the Planck mass; I -C Mp, as is true in most infiation models, and (iii) the field is also 
much smaller than Mp; (j) -C Mp. The last two conditions allow us to make a perturbative 
expansion of the potential in terms of / and cj). Under these assumptions, we write the most 
general Kahler potential and superpotential : 



K 



W 



Ml 



I^I(f)^(f) + 



d 
Ml 



I^IW{(j)) + h.c. 



(10) 



neglecting terms with higher powers of I^I, W{(f)) which give more suppressed contribu- 

tions. The conditions (ii) and (iii) imply that Fj ^ diW ^ V^HMp and W/Mp = alH < Fj 
with a < 1. Then, we have 



;i-c), ^ V3 [l - c(n + 3) - V3i 



Mr 



H 

Mp 



(11) 



As advertised, one obtains \C^\ ~ 1 but \Ca\ <^ 1 |H1- 

If H replaces m</,, Eq. Q becomes v ~ (ifMp)^/""*"^, which leads to the following; 

f (i/M^)i/"+^ during infiation, 

f~< [max(if, m3/2)Mp]^/"+^ after infiation, 
[ (^«3/2^p)^^"^^ = fpQ after reheating. 



(12) 
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While the radial degree of freedom follows such a thermal history, the corresponding phase 
degree of freedom, a, will enjoy the following potential of the pseudo Nambu-Goldstone boson 
(PNGB) type 0: 

(13) 









\Mp) 


1 — cos 









V{a) ^ {CaH + Ay 
which is very flat during and after inflation. From Eq. ()12|). one finds 

{CaH^ during inflation, 

max(C^iJ^, AiJ, Am3/2) after inflation, (14) 
Ams/2 after reheating, 

where we assumed that, during inflation iJ* > A/Ca- 

4 The curvaton cosmology 

It is now obvious that the field a will naturally realize the curvaton paradigm. Taking Eq.(jH} 
literally, a can be parameterized as 6* = a /\/2v, that is, (j) = ve"^ [cf. Eq. (jHl)]. To obtain simple 
predictions, let us assume that the value of a is completely randomized during inflation, so that 
in our part of the Universe 9 at horizon exit has some value < 27r. (We will come back to 
this point later). 

It remains constant as far as H ^ rricr and the oscillation starts at if ~ 7713/2- Finally, the 
curvaton decays and thermalizes at if ~ Fg- which is determined to be quite late due to its 
axion-like couplings. The interaction of a with ordinary particles is governed by the effective 
/i term in Eq. ^ leading to /i = hfp^/Mp ~ m3/2. Thus, the decay rate of a into two Higgs 
field is 



where we put /x = rria = 7723/2. Taking the approximation of ~ ^3/2/ fpQ fpQ ~ 
(7?73/2Mp)-^/""'"^, we get the decay temperature of a as follows; 

r 160 GeV for 7i = 1 
Tdec ~ I 0.6 GeV for 71 = 2 (16) 
I 30 MeV for 71 = 3 

with 777.3/2 = 300 GeV. Here we have restricted ourselves to < 3 for which the decay occurs 
well before the nucleosynthesis occurring around Tbbn ~ 1 MeV. 
The density perturbation driven by the curvaton is given by [2] 

C-^^^ (17) 
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where r is the density ratio of p^/ p-y at the time of the a decay: 

1 rR^V/' r a, \2 



6 V r„ y \M, 



;i8) 



Here Hm = min(mo-, F/) where F/ is the inflaton decay rate determining the reheating tem- 
perature. Note that, in our model, the values of r, H^, and a* are determined by dimensionality 
of the superpotential in Eq. (0), which are related to the axion scale, fpq. 

Inserting Eqs. ()12j] and (fTH|l to Eqs. p7|l and (fTK|) . we can determine the inflation Hubble 
parameter H^, and the density ratio r which reproduce the observed value of C ~ 5 x 10^^. 
First, the value of r can be calculated from the relation: 



r V AnC e^+' ( H^Y ( Mp\ ^("+1) 



l + |r 6"/2 V"^3/2/ \"^3/2/ 

and then it fixes the inflation Hubble parameter as 

n+l 

|..,4.c«,,=*^(i±l:)". (.0) 

For a typical set of input parameters; ^3/2 = 300 GeV and Hm = 1 GeV, we get the following 
values for H^/Mp and r; 



( 



X 10-^^0"^/^ , 1.4^0^^) for n = 1 

, r) ~ <! (1 X 10-5 ^0^^ , 6 X 10^ e^) for n = 2 (21) 



(4 X 10^5 ^0''^ , 2 X 10^ ^0^^) for n = 3 



which is valid for the typical case of 60 < 1, i.e. r < 1 for n = 1 and r ^ 1 for n = 2, 3. The 
exact behaviors of {H^../Mp,r) depending on 6q are shown in Figures 1, 2 and 3. 

The value of r has to be large enough to avoid excessive non-Gaussianity in the density 
perturbation spectrum. The constraint from the WMAP observations reads: r > 9 x 10~^ (23- 
Enforcing this constraint in Eq. (fT^ we obtain the bound 

0.232 f H^Y^ f MpY''^^'^' 
(103.37rC)V("+i) 1,^3/2^ 1^3/2^ ■ ^ ' 

Suppose first that we live in a typical part of the Universe, then 6q ~ 0{1) (see appendix 
for a discussion on the randomization of the curvaton field). In this case the value of r is always 
large enough to avoid excessive non-Gaussianity. For n = 2, 3, H^, becomes barely compatible 
with the limit H^/Mp < 10~^ coming from the observational bound on the primordial tensor 
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perturbation. (Actually, the curvaton model requires a somewhat smaller limit [24J, assuming 
that inflation is of the slow-roll type.) The case with n = 1 then gives most satisfactory result 
with ~ 10^^ GeV. In this case, the curvaton decay temperature is T^^c ~ 100 GeV and the 
resulting entropy dumping is not significant (it is comparable to the preexisting one), so that 
the usual cosmological properties are retained. That is, the conventional dark matter candidate, 
the lightest neutralino or the axion with fpq = 3 x 10^° GeV, is equally acceptable within our 
curvaton context, and the usual baryogenesis mechanism can also work. Note, also, that in this 
case we can have substantial non-Gaussianity in the density perturbation spectrum. Indeed, 
from Eq. PT|) we readily obtain j2I] 

/nl = ^^0.9^o"'^^ (23) 

The largest value of the above corresponds to the smallest possible value of which is deter- 
mined by Eq. (j2DI), when we demand that H^/Mp < 10~^. With a little algebra it can be shown 
that (^^o)mm — 0.34:y/r, which gives (/NL)max — O{10). Such a value of /nl may be observable 
in the near future. 

In the above we assumed a reheating temperature Treh ~ lO^GeV, which saturates the con- 
straint of gravitino overproduction. This constraint is important, in the case when n = 1, 
because the entropy production by the curvaton decay is not enough to dilute the gravitinos. 
However, a lower Treh corresponds to a smaller Hm, which results in a smaller r, according to 
Eq. p9|) . This, in turn, results in a larger [cf. Eq. ()20|)]. which again violates the upper 
bound on H^,. On the other hand, a higher Treh corresponds to the range 1 GeV < Hfh < m3/2, 
which, for a given n, can increase r somewhat. If r ^ 1 then the entropy production by the 
decay of the curvaton relaxes the gravitino constraint. However, as shown by Eq. (j^Uj) . loses 
its sensitivity to a large r, even though, as can be seen in Figure 1, may be reduced down 
to 10^'^Mp for Hffi ~ m3/2.^ 

Alternatively, we may allow the possibility that 6q is untypically small in our region. Then 
the value of if* needed to generate the observed density perturbation is reduced by a factor 

n+l 

6q " for n = 2, 3. Now the cases n = 2, 3 become viable, but the n = 1 is not because it violates 
the if* < lO^^Mp bound [cf. Eq. Moroever, for n = 2 and 3 we now have the fascinating 

possibility that the curvature perturbation has the observed magnitude only because we live 

■^Note, here, that there is also a lower bound on _ff*, as was derived in 125 . In our case it can be shown 
that, because v during inflation is different that in the vacuum, this bound reads > {v^/ fpQ)*^'^10-'^^Mp 
|2Sj. Using, Eqs. © and ^ this bound can be recast as iJ* > [IQ-'^^^^+'^^Mp /m3/2)]^ Mp. The tightest 
case corresponds to n = 1, which gives > 10~^Mp. Clearly, our results in Eq. (|^ satisfy well this bound 
for 6*0 - C(l). 
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in a special part of the Universe. (To be precise, the bound < 10~^Mp makes this state of 
affairs inevitable for n = 2, 3) This is actually a quite common outcome of curvaton models jH], 
and provides a good example of how a theoretical model can make anthropic considerations 
mandatory (in absence of a model to fix 6 dynamically to a small value). 

In this model, Peccei-Quinn symmetry is broken during inflation and is not restored. No 
axionic strings are formed, and axions are produced by the oscillation of the homogeneous axion 
fleld. Their density is jTH] 




= (1 to 10) — —' D, (24) 



where N is the multiplicity of the axion vacuum and D is the dilution coming from entropy 
production after axion creation. (In our model, D = 1 for the most interesting case n = 1.) 
The uncertainties are big, but it is clear that the axions can be the dark matter (fia = 0.25) 
for roughly of order 1. 

The perturbation will generate a matter isocurvature perturbation S^, which should be 
< 0.1( to be compatible with observation. It is given by 

Sm = (25) 

3 Pa 

But 66a ~ 600 because both are generated from the vacuum fluctuation. Therefore, in view of 
Eq. (Hlj), we flnd 

which can be < 0.1 as required by present observation. On the other hand it could be observ- 
able in the future. 



5 Discussion and Conclusions 

Within the curvaton paradigm, we have proposed a new candidate for the fleld which causes the 
curvature perturbation. We believe that it is one of the most attractive candidates yet proposed 
for that fleld. Being a fleld that is part of the flaton realization of PQ symmetry-breaking, and it 
can easily be kept light during and after inflation. In a typical part of the universe, the observed 
magnitude ~ 10~^ of the curvature perturbation comes from a modest and quite reasonable 
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hierarchy between the inflationary Hubble parameter and the vacuum expectation value of the 
Peccei-Quinn field. The predicted values of the axion scale and the Hubble parameters are 10^° 
and 10^^ GeV, respectively, corresponding to = 1. In this case, the entropy dumping due to 
the curvaton decay is negligible so that the conventional cosmological predictions concerning 
dark matter components and baryogenesis remain unchanged. The model can generate two 
kinds of isocurvature pertubations. One is an axion isocurvature perturbation, uncorrelated 
with the curvature perturbation. The other is an isocurvature perturbation in some other 
kind of dark matter, or in the baryonic matter, which is produced either before or at curvaton 
decay. As described in j^H I2H1, such a perturbation is generic to curvaton models, and is 
fully correlated or anti-correlated with the curvature perturbation.^ The detection of any such 
isocurvature perturbation would be evidence in support of the curvaton model. 

Finally, we have shown that, in this case, non-Gaussianity at at an observable level (with 
/nl ~ C^(IO)) is possible. Our curvaton candidate is a PNGB, corresponding to an angular 
degree of freedom associated with the QCD axion. As a consequence, it can avoid the usual 
mass of order H during (and after) inflation, because supersymmetry breaking affects the mass 
only through A terms, which can be controlled by appropriate symmetries. 

Our model of PQ symmetry breaking is similar to one already extensively investigated 
|2T| l29] . except that we take the PQ symmetry to be spontaneously broken throughout the 
history of the Universe. In contrast, the investigations of pT| EH] assume that the symmetry 
is initially unbroken, leading to thermal inflation. These possibilities correspond respectively 
to radial masses- squared of order ^H"^, and the two signs should be deemed equally likely in 
the absence so far of a string-theoretic prediction. The unbroken paradigm has a very different 
cosmology from the one that we are adopting, producing in particular copious saxion- or axino- 
like particles which may decay into relativistic axions whose energy density is enough to affect 
nucleosynthesis • Also, the lightest axino-like particle may be the lightest supersymmetric 
particle pS], providing a more natural implementation of the Cold Dark Matter scenario that 
was originally proposed [20] in the context of non-flaton models. None of this occurs within 
our paradigm. 

As it invokes the axion, our model, as it stands, is open to the criticism that the axion 

mass is implausibly small. Indeed, in the kind of models that we have considered where the 

axion is an angular part of a complex field, a non-renormalizable term in the potential with 

dimension d generically breaks PQ symmetry and contributes to the axion mass an amount 

Correlation or anti-correlation of the CDM {baryon} isocurvature perturbation depends on whether CDM 
creation {baryogenesis} takes place (just) before curvaton decay or due to the curvaton decay itself. Note, that 
there is no residual isocurvature perturbation if CDM creation {baryogenesis} occurs after curvaton decay. 
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of order ~ y^<^-'^)l'^ in Planck units. To keep the axion mass to the required value of order 
10~'^°Mp, terms up to d ~ 12 must respect the PQ symmetry. (We take v ~ fpq = 10^^ GeV 
for an estimate.) A widely-discussed possibility to ensure this is to impose a Zn subgroup of 
the PQ symmetry [32].^ 

Because of these considerations, one may favour a solution of the CP problem in which the 
QCD axion is identified with a string axion [HI], or one without any axion at all.^ In that case 
the field considered in our proposal becomes ad hoc, introduced solely to explain the curvature 
perturbation. Still, because of its simplicity and the ease with which it is made sufficiently light, 
we feel that the present curvaton model is extremely attractive in comparison with others. Let 
us end by mentioning those other models which have a close connection with ours. 

The closely related models are considered in Refs. jElEI, where as in our case the effect 
of the angular part of a flat direction has been considered. The difference from our case is 
that the flat direction is supposed to have zero VEV, so that the curvaton does not exist as a 
particle in the vacuum. As a consequence, the phase can only induce the curvature fluctuations 
in the radial field which leads to different predictions. In addition, the proposal of jT2I does 
not invoke the mass-squared of order —H^, so that the radial potential is very flat with only 
a two-loop thermal correction breaking the symmetry. (In both cases, successful curvature 
perturbations can arise for n = 3, contrary to our case with n = 1.) Because of such features 
the success of these models depends on computations which are much more tricky than in our 
case, though that is of course not necessarily an argument against them. However, we note 
that the computation of ^21 5 involving the very flat radial potential, works only if inflation 
lasts for a limited amount of time, because it takes the radial field during infiation to be at the 
edge of the slow-roll regime, V" ~ H^. It seems more reasonable to assume that infiation lasts 
long enough to allow the quantum fiuctuation to randomize the radial field within the smaller 
region V < H^. (The analogous assumption for the angular field is of course the one that we 
made.) 

^Imposing the full PQ symmetry of course works, but exact continuous global symmetries are widely regarded 
as incompatible with string theory and even the existence of gravity. A different possibility, which does not 
seem to have been mentioned before, might be to suppose that the PQ fields are actually moduli, making 
the origin a point of enhanced symmetry. Then the non-renormalizable terms may be suppressed by a factor 
(TeV/Mp)2 ^ 10^3°. However this is small enough for all 0? > 5 only if fpq has an implausibly small value of 
order 10* GeV. In any case, this mechanism cannot be used in the flaton case, which invokes an unsuppressed 
coefficient for one term. 

^The string axion as a curvaton candidate is discussed in 0. As the string axion is a PNGB one may 
be hopeful that it can be kept sufficiently light in the early Universe, but string phenomenology is not yet 
sufficiently developed that one can be sure. 
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The other related model ^I] also invokes a fiaton model of PQ symmetry breaking, but 
now the curvaton candidate is a radial field and it is not clear how to keep it light in the early 
Universe. Also, the model works only if infiation is of short enough duration that the curvaton 
does not enter the randomization regime. 

Acknowledgments: EJC was supported by the Korea Research Foundation Grant, KRF- 
2002-070-C00022. 

A Appendix: Randomization of the PNGB curvaton 

In this appendix we elaborate more on the assumption of the complete randomization of the 
value of the curvaton during infiation. Since the curvaton is an effectively massless field dur- 
ing infiation, its value is perturbed by the action of quantum fluctuations, which introduce a 
perturbation of the order of 5a = Th per Hubble time, where Th = H^/27c is the Hawking tem- 
perature of de-Sitter space. Provided that the potential V{a) is not too steep, these quantum 
'kicks' can move the fleld around by random walk so that, given enough e-foldings, the phase 
6q may assume an arbitrary value (typically of order unity) by the time when the cosmological 
scales exit the horizon. The criterion for this to occur for a PNGB curvaton is > He, where 



The above critical value has been obtained by considering the fact that the region of V{a), 
where the fleld is randomized, is determined by the condition V{a) < This condition 

can be understood as follows: The kinetic energy density of the quantum 'kicks' is ~ H^. 
Hence, the quantum fluctuations can displace the fleld from its minimum only up to potential 
density V ~ T^. For a PNGB V{a) ~ m'^a'^, which means that the borders of the randomized 
region are at aq ~ / . If if* is large enough then the randomized region includes the entire 
range of a. Setting aq ^ v we obtain the critical value of H^, shown in Eq. (j28|) . over which 
the PNGB is fully randomized. 

Using Eqs. (fT^ and (fT^ the bound ii* > He translates into 



where we also used Eq. (1201) with r > 1, according to Eq. (|2ip. With 6'o ~ we see that 

complete randomization of cr* requires Ca < 10"''. According to Eq. (fTT|) . this condition is 
satisfled provided that |/| < 10~^Mp. ^ 

^Note that, in this case, there is no problem with Ca > A/H^ ~ 10^^°. 



He = y/rUaV. 



(28) 



2n 




(29) 



12 



If complete randomization is not realized then a^, may be much smaller than v during 
inflation. In this case, according to Eq. (fT^ . r can be smaller than unity and the results in 
Eq. (j?H) are modified. In particular, the typical value of is estimated as 

H;/m^ for < He 

where we considered that, in the non-randomized case, the typical value is a* ~ ag. Inserting 
the above into Eq. p7|) and after some algebra it is easy to find 



r nr HJH^ for H^> 

4 + 3r ^Ym, [H,/H, for < H, 

It it easy to see that f(r) is an increasing function of r. Hence, the WMAP non-Gaussianity 
constraint reads: f(r) > f(9 x 10~^) ~ 2 x 10"'^. For a successful PNGB curvaton, this con- 
straint has to be satisfied at least for the maximum allowed f. Therefore, in view of Eq. ()3ip . 
we obtain 

rmax ^ <^v/m„ > 2 X 10^1 (32) 
Using Eqs. ((12)), and ((201), we get 

'"max — rrf~ a ^ "^ma.^ < n ' V"^"^/ 

Combining Eqs. (jH^ and (jH^ we obtain 

which is typically satisfied by a successful curvaton. 
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Figure 1: Log-log plots of H^/Mp with respect to r for all cases: n = 1,2,3. The dashed line 
corresponds to the bound: H^/Mp < 10^^. The upper graph corresponds to 6*0 ~ 1, whereas 
the lower graph corresponds to ~ 0.1. We see that, in the former case, only the case of n = 1 
manages to avoid the bound. On the other hand, in the latter case, all three possibilities are, 
in principle, allowed. The curves meet when H^, ~ Mp as expected. 
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^0 



Figure 2: Log- log plot of 60 with respect to r for all cases n = 1, 2, 3. The solid lines correspond 
to Hm ~ 1 GeV, whereas the dotted lines correspond to Hm ~ ?tJ'3/2- The horizontal dashed 
line corresponds to ^0 ~ 0.1. We see that, for > 1 GeV, the value of r is slightly larger, 
for a given 6*0 (it can increase up to an order of magnitude at most). Obviously, 6*0 is bounded 
from above as log6'o < logvr ^ 0.5. 
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Figure 3: Log- log plots of H^/Mp with respect to 9q for all cases: n = 1, 2, 3 with Hf^ ~ 1 GeV. 
The dashed line corresponds to the bound: H^^/Mp < 10^^ (which ensures the compatibility 
with the observations), while the dotted line corresponds to the bound: H^/Mp < 10^^ (which 
ensures that the inflaton does not contribute significantly to the curvature perturbations). We 
see that, when Oq ~ 1, only the case of n = 1 is possible with H^/Mp ~ 10~^, which means that 
the inflaton generated curvature perturbations are not negligible. For the curvaton to be the 
sole source of the curvature perturbations one has to limit oneself to the n = 2, 3 cases, where 
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